Valley- dependent Brewster angles and Goos-Hanchen effect in strained graphene 
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We demonstrate theoretically how local strains in graphene can be tailored to generate a valley 
polarized current. By suitable engineering of local strain profiles, we find that electrons in opposite 
valleys (K or K') show different Brewster-like angles and Goos-Hanchen shifts, exhibiting a close 
analogy with light propagating behavior. In a strain-induced waveguide, electrons in K and K' 
valleys have different group velocities, which can be used to construct a valley filter in graphene 
without the need for any external fields. 
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Introduction. In recent years, graphene, a single layer 
of carbon atoms arranged in a hexagonal lattice has 
shown abundant new physics and potential applications 
in carbon-based nanoelectronic devices. p]-[3| The novel 
properties arise from the linear energy dispersion and the 
chiral nature of electrons at the K and K valleys of the 
Brillouin zone. [3| Mathematically, the Dirac equation 
describing the motion of massless quasiparticles is very 
similar to the Helmholtz equation for an electromagnetic 
wave, except for the negative energy spectrum. This re- 
markable fact could make it possible to observe optical- 
like electron propagating behavior and to construct Dirac 
electron-optical devices with graphene. [4] The mean- free 
path of electrons in graphene can approach microns at 
room temperature, making electrons behave ballistically 
in a graphene microstructure. [3, Q In such a ballistic 
regime, the scattering of electrons by potential barriers 
can be understood by comparing with the reflection, re- 
fraction and transmission of electromagnetic waves in in- 
homogeneous media. Therefore, graphene is a testbed 
to examine optical-like phenomena of Dirac fermions. 
Several studies have been devoted to such optical-like 
electron behaviors using external electric bias. 0, 04HI 
These proposals require to deposit a metallic gate above 
graphene, which is an additional complication. 

The strain effect in graphene provides a new way to 
manipulate electron transport without external fields. 
The two valleys in graphene are related by time reversal 
symmetry and act in much the same way as electron spin 
in spintronics. The valley degree of freedom in graphene 
might be used to control the characteristics of graphene- 
based devices, referred to as valley tronics. [16l-ll8| Re- 
alization of such valleytronics requires an effective and 
robust scheme to produce valley polarized currents, al- 
though several schemes about valley filters were proposed 
utilizing, e.g., the edge profile of graphene nanoribbons, 
the trigonal warping effect in the energy spectrum, and 
intense irradiation on a bilayer graphene. [18l-[2H| In this 
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letter, we propose a simple and effective way to realize 
some optical-like behaviors of electrons and valley filters 
by using the strain effect alone. We find that electrons in 
opposite valleys can be perfectly transmitted or totally 
reflected in the presence of strain. A quantum waveg- 
uide formed between two strained regions (see Fig.[T]) can 
confine electrons in it and the reflected beam is shifted 
laterally along the interface by a distance with respect to 
the incident beam, as the Goos-Hanchen (GH) effect [9] 
in optics. We demonstrate that the valley-dependent GH 
effect in graphene results in different group velocities for 
electrons in K and K' valleys. 

Model. Mechanical strain in graphene can be described 
by a gauge vector potential. To illustrate the operat- 
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FIG. 1: Schematic of a strained four-terminal guiding de- 
vice including contacts I (injector), C (collector), and Gl, 
G2 (electrical ground 1, 2). Strain exists in the shaded re- 
gions. The black arrow denotes the unpolarized incident elec- 
tron beam. The red (blue) arrows stand for the K- valley 
(if Valley) polarized electron beams. A valley detector as pro- 
posed in Ref. [22j] could be connected to contact C to measure 
the polarization. The detector is also a graphene-based device 
with a substrate strain and a ferromagnetic (FM) strip. 
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FIG. 2: The transmission probability T as a function of the 
incident angle <j>. In the strained region, the wave vectors of 
electrons in K (K 1 ) valley satisfy k 2 x ^ + (k y + t;A y ) 2 = E 2 as 
indicated by the red (blue) arc curve. The red (blue) dashed 
region indicates the transmission window for electrons in the 
K (K f ) valley, (a) E = 10, A y = 5, D = 1. (b) E = 10, 
A y = 15, D = 2. 



ing principle of the proposed device, we take a simplified 
gauge field which is equivalent to a delta-function-like an- 
tiparallel magnetic double barrier perpendicular to the 
graphene monolayer, but with opposite direction at K 
and K' valleys, [Il,^ The Landau gauge A = (0, A y ,0) 
is adopted in our calculation. The length scale of the spa- 
tial variation of the pseudo gauge vector potential is much 
larger than the lattice spacing of graphene, which implies 
that inter valley scattering is weak at low-energy regions. 
The low-energy electrons can be well described by the ef- 
fective Hamiltonian H = v F a^ ■ (p+^/vp) + V\ 
where the superscript i indicates the different regions (see 
Fig. [2]), vf is the Fermi velocity, cr^ (a' = —a ) are 
the pseudospin Pauli matrices, p is the electron momen- 
tum, V 1 is the electrostatic potential in the region i, and 
£ = =pl labels K and K' valleys. In the calculation, 
we introduce the dimensionless length and energy units: 
Ib = (h/eBo) 1 / 2 , Eo = Tivf/Ib (which are 81.1 nm and 
7.0 meV for a typical pseudomagnetic field Bo =0.1 T). 

Brewster-like Angle. We start by investigating electron 
transmission through a region of uniform uniaxial strain 
with width D as shown in the inset of Fig. [2j Here, we 
set V 1 = for simplicity since the key physical mecha- 
nism to introduce the valley-dependent Brewster angle is 
the gauge vector £A. The translational invar iance along 
the y direction gives rise to conservation of transverse 
wave vector k y , and thus the solutions can be written 
as ifj(x,y) = ip(x)e' lkyy . In the strained region, the lon- 
gitudinal wave vector k x ^ = (E 2 — (k y + ^A y ) 2 has a 
different dependence on the vector potential for the two 
valleys denoted by £ = =pl, since the strain-induced pseu- 
domagnetic fields have different signs for the two valleys. 
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where the incident and refractive angles are defined as 
cj) = arc sin (k y /E) and 0% = ar c sin [(k y +£A y )/E]. These 
expressions are valid for incident angle <\> < (j) c ^ = 
arcsin(|A: y + £A y \/E). <p c ^ is the critical angle for total 
reflection in optics and sin 6^/ sin^> = (k y +t;A y )/k y = n 
gives Snell's law for transmitted electrons. Note that 
when k y (k y + £A y ) < 0, the refractive index n of the 
strained graphene is negative just like for a metamate- 
rial with a negative refractive index. Importantly, one 
can tune the refractive index n mechanically in quite 
a large range, which is not so for the metamaterials. 
When <j) > C? £, electron beams can be totally reflected at 
such a strained barrier, since the longitudinal wave vector 
k X £ becomes imaginary accounting for the occurrence of 
evanescent modes in the barrier. 

Figure [2] shows the angular dependence of the trans- 
mission probability = 1 — |r^| 2 , which exhibits a re- 
markable valley-dependence. It is obvious that T^(0) = 
TJf(— </>), where £ = — £, i.e., the transmission of electrons 
in K and K r valleys shows a mirror symmetry about 
(j) = 0. This feature is guaranteed by the time reversal 
symmetry in strained graphene. In the presence of uni- 
form uniaxial strains, the wave vectors of electrons in K 
valley satisfy the relationship k 2 _ x + (k y — Ay) 2 — E 2 
as indicated by the red arc curve. When \k y — A y \ > E, 
(equivalent to G [— 7r/2, c ,-]), the wave vector k x in 
the strained region becomes imaginary, indicating the ap- 
pearance of evanescent modes, thus transmission is to- 
tally blocked. The critical angle for total reflection is 
given by Snell's law, <p c - = arcsin^A^ — E)/E). Thus 
the transmission window of electrons in K valley is re- 
stricted to the red shadowed region (cj) G [0 Cj _, 7r/2]), 
as shown in Fig. [2](a). In contrast, the transmission of 
electrons in K' valley shows mirror symmetry behavior 
with respect to that in K valley. Similarly, the criti- 
cal angle for electrons in K' valley is given by </> C5+ = 
arcsin((£^ — A y )/E). As a consequence, there exist mir- 
ror symmetric transmission windows for K and K' elec- 
trons that can be completely separated by increasing the 
strain gauge fields [Fig. [51(b)]. One can see in Fig. [2] that 
for some specific incident angles 0b, resonant peaks that 
reach perfect transmission exist in different windows for 
electrons in K or K' valley. This can be readily checked 
by analyzing the zeros of the reflected amplitude r corre- 
sponding to sm(k x ^D) = 0. The perfect tunneling peaks 
are fully coincident for k x ^ = nn/D [see Fig. [21(b)]. The 
number of resonant peaks increases with the incident en- 
ergy E and the width D of the strained region. Note 
that if such resonant peaks for one valley £ are located in 
the transmission gap for the other valley £, we can obtain 
valley-^ polarized transmitted electron beams and valley- 
£ polarized reflected beams. These characteristic angles 
are analogous to the so-called Brewster angle in optics, at 
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FIG. 3: The energy spectra for the lowest channel modes in 
the stain-induced waveguides for (a) A y = 5, D — 1, and (b) 
A y = 15, D = 2. 



which s- and p-polarized light are produced. We can pro- 
pose a valley polarization device that is straightforward 
and analogous to the Glan-Thompson optical polarizers. 
The strained stripe leads to perfect transmission of one 
valley component and total reflection of the other val- 
ley components when electron beams are incident at the 
Brewster- like angles <\>b> 

Valley- dependent Goos-Hanchen Effect. In analogy 
with total reflection of light at the interface between 
two materials with different refractive indexes, the elec- 
tronic GH effect describes the shift gqh of the reflected 
electron beam for total internal reflection at an inter- 
face, along the transverse direction. The GH shift could 
be positive or negative which is determined by external 
electric [To[ and/or inhomogeneous magnetic fields 
an analogue with the optical GH effect in normal ma- 
terials or metamaterials. A uniform uniaxial strain 
in graphene can lead to a valley-dependent GH shift 
o gh,£ — (2k y + £A y )/(k x K£) : where k is the modulus of 
the imaginary wave vector in the strained region in the 
case of total internal reflection. For a waveguide formed 
in between two strained stripes as shown in Fig. [TJ the 
transmission gap at each interface can lead to the confine- 
ment of electrons in the channel associated with multiple 
total internal reflections. The difference between the tun- 
neling forbidden regions for electrons in different valleys 
gives rise to valley-dependent electron propagation along 
the channel. In the channel, the electronic GH shifts at 
each interface are accumulated during multiple reflection 
processes. 

The GH effect is found when there is total internal re- 
flection, and thus there are bound states localized in the 
channel between the two interfaces as shown in Figs.[3ja) 
and [3Kb). The bound states exist only in the region be- 
tween the curves given by the equations \k y + £^1 = \E\ 
and \k y \ = \E\, i.e., the transmission gap as we discussed 
above. The spectrum of each valley (K or K') alone 
is asymmetric, but for both valleys (K and K'), the 
total spectra are mirror symmetric with respect to the 
transverse wave vector k y , as required by time reversal 
symmetry. That is very different from bound states in 
a pure electric waveguide in graphene, where the bound 
states are valley-independent and always mirror symmet- 



ric with respect to k y = 0. [To| The number of bound 
states increases with the strain-induced gauge filed A y . 
Three typical regions in the energy spectrum are found: 

i) E < Ay/2, i.e., (j) C}T = ±7r/2, the incident electrons 
from both valleys are always totally reflected at the in- 
terfaces. The bound states for K- or K'- valley electrons 
coexist for all incident angles. The electrons in differ- 
ent valleys have different group velocities determined by 
the slope of the energy dispersion relation. Thus such a 
waveguide may be used to separate the electrons in dif- 
ferent valleys after passing a sufficient long channel. No- 
tice that there are several local minima in the dispersion 
relation indicating a vanishing group velocity. Thus elec- 
trons in K or K' valley will be trapped in the waveguide, 
which can be used to construct a valley memory device. 

ii) Ay/2 < E < A y , i.e., ir/2 > Cj _ > (/> c , + > -7r/2, the 
incident electrons from both valleys are confined in the 
channel for incident angles <\> G [0 Cj +,0 Cj _]. The trans- 
mission in such an angular window is similar to that in 
situation i). For <p < C?+ (or <j> > C ,-), only the bound 
states for K valley (or K' valley) electrons appear, while 
electrons in K' valley (or K valley) are able to penetrate 
into the strained region and eventually disappear from 
the channel region. So we can produce a valley-polarized 
current in such a strain-induced waveguide, iii) A y < E, 
i.e., 7r/2 > C?+ > Cj _ > — 7r/2. In this case, there 
is no coexistence of bound states for K or K' valley. 
The transmission properties are similar as that in case 
ii) for <j) < (/> c?+ and <\> > 4> c ,-- Note that this valley- 
dependent transport property can not be realized in nor- 
mal electric/magnetic waveguides proposed in previous 
works (see Refs. [13, El). 

Finally, we focus on how the valley-polarized transmis- 
sion features shown above are reflected in the conduc- 
tance G. Recently, Williams et al. [15| reported fiber- 
optic guiding of electrons in graphene. Differential re- 
sistance can be measured using a standard lock-in tech- 
nique and thus guiding efficiency can be extracted. Our 
proposal depicted in Fig. [1] is a valley-related version 
of their setup. Here the guiding modes serve as valley- 
filtered modes while the bulk modes are unpolarized. To 
distinguish these valley-filtered modes from bulk modes, 
a higher guiding efficiency is desired. High guiding ef- 
ficiency can be achieved by applying a uniform mag- 
netic field as reported in Ref. [15]. But in this letter we 
propose a different scheme, where an applied mechani- 
cal strain behaves like a delta-function-like pseudomag- 
netic field, with opposite direction at the K and K' val- 
leys. We performed numerical simulations of electrical 
conduction in a tight-binding model [19] of a strained 
graphene stripe. The dispersion relation E(k y ) is plotted 
in Figs. [U a) and 2Kb). The guiding modes are formed 
in the interval of the bulk states as shown in the dashed 
pentagon in Fig. 2Kb). Note that when E < 3 eV, all the 
forward-going guiding modes are located in K' valley, 
while the forwarding-going modes in K valley are solely 
bulk modes. The density distributions of electrons in the 
lowest two forward-going modes in K and K' valleys are 
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FIG. 4: (a) Energy dispersion of a strained graphene rib- 
bon with width N=140 (in unit of honeycombs), pjj] Lateral 
strains are applied on the ribbon except the middle part as 
shown in Fig. [T] Ko is a reciprocal lattice vector, (b) A 
blowup of the box region in panel (a), (c) Density distribu- 
tions of electrons in the lowest two forward-going modes with 
k y = — O.GOifo in K' valley and k y = — 0.18i^o in K valley, 
(d) Conductance versus incident energy. The blue (green) 
curve gives the number of propagating modes in the channel 
(lead). 

shown in Fig. ffl^c). Electron guiding modes (in K' val- 
ley) and bulk modes (in K valley) are separated spatially 
in the channel and the strained region, respectively. The 
conductance shows good correspondence with the num- 
ber of guiding modes rather than the much larger number 
of transmission modes in the leads [see Fig.HJd)]. This is 
because electrons in bulk modes in the lead leak out of the 



channel and eventually disappear together with those in 
bulk modes in the strained region via the electric grounds 
Gl and G2. The reason for the oscillation of the conduc- 
tance is that the stripe has a finite width and some bulk 
modes in the strained region still affect the conductance. 
One would expect that the conductance will approach 
the guiding mode steps very well when the strained re- 
gion are electrical grounded. The guiding modes can ef- 
fectively carry current from the injector I to the collector 
C and thus dominate the channel conductance. As we 
demonstrated that the guiding modes are valley- filtered, 
the desired valley polarized current can be measured in 
collector C. A valley detector as proposed in Ref. [22| can 
be connected to the contact C to identify the polariza- 
tion. The detector shows a large valley-resistance ratio, 
in analogy to the giant magnetoresistance effect. Thus 
valley polarized current could be identified by measuring 
the conductance at terminal C2. 

In summary, We demonstrated theoretically how a 
strain-induced gauge field can be tailored to generate 
valley-polarized transport in a single layer of graphene. 
Our results show that the propagating behavior of elec- 
trons in graphene exhibits deep analogies with light as 
a consequence of the massless linear Dirac dispersion 
and the chirality of electron states. Most importantly, 
the gauge fields induced by strain can lead to valley- 
dependent transport phenomena, e.g., the Brewster an- 
gles and the Goos-Hanchen effect. Electrons in a waveg- 
uide structure formed by two strained stripes propagate 
along the channel with different velocities at different 
valleys. This feature shed new light on constructing 
graphene-based valleytronic device by using only strain, 
without the need for any external fields. 
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